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In this work, we study the isospin-violating decay of <j> — * urn and quantify the electromagnetic 
(EM) transitions and intermediate meson exchanges as two major sources of the decay mechanisms. 
In the EM decays, the present datum status allows a good constraint on the EM decay form factor 
in the vector meson dominance (VMD) model, and it turns out that the EM transition can only 
account for about 1/4 ~ 1/3 of the branching ratio for <f> — > loty . The intermediate meson ex- 
changes, KK(K*) (intermediate KK interaction via K* exchanges), KK*(K) (intermediate KK* 
rescattering via kaon exchanges), and KK*(K*) (intermediate KK* rescattering via K* exchanges), 
which evade the naive Okubo-Zweig-Iizuka (OZI) rule, serve as another important contribution to 
the isospin violations. They are evaluated with effective Lagrangians where explicit constraints 
from experiment can be applied. Combining these three contributions, we obtain results in good 
agreement with the experimental data. This approach is also extended to J/ip(ip') — > uir°, where 
we find contributions from the KK(K*), KK*(K) and KK*(K*) loops are negligibly small, and 
the isospin violation is likely to be dominated by the EM transition. 

PACS numbers: 12.40.Vv, 13.20.Gd, 13.25.-k 



I. INTRODUCTION 

The isospin breaking decay channel 4> — > u>ir° has been measured by experiment with improved pre- 
cisions 1], and the Particle Data Group quote BR{<j) — > urn ) — (5.2~t_\'i) x 10~ 5 as the world average 
for its branching ratio @|. This decay channel is very interesting due to the presence of the OZI-rule 
violation and isospin symmetry breaking together. These two mechanisms, which generally account for 
different aspects of the underlying dynamics, are correlated in this channel. With the available of much 
improved experimental information about other related transitions, one can pursue a quantitative study 
of the underlying dynamics and learn more about the correlation between the OZI-rule violation and 
isospin symmetry breaking in the non-perturbative regime. 

The electromagnetic (EM) decay of 4> — > ojtt is an important source of isospin violations, where the s 
and s annihilate into a virtual photon, which then decays into ojtt . The other source of isospin violation 
originates from the mass differences between the u and d quark Q. It can contribute to <j> — > ujit via 
OZI-rule- violating strong decays. 

In the literature the isospin violation in 4> — > ojtt was studied by isoscalar and isovector mixing, e.g. 
<f>-oj-p° and rj'-rj-ir mixings [4|, la [y, [7|, 13, l2( . This scenario contains both EM and strong transitions in an 
s-channel, and allow the <j> — » ujit decay without violating the OZI-rule [l(l EH ■ In such an approach, the 
EM and strong decays cannot be separated out. An alternative view is to separate the EM and strong 
processes by explicitly introducing the EM amplitude as an s-channel process, and then including the 
hadronic loop contributions as the f-channel processes. This will be our focus in this work. Our strategy 
is to constrain the EM transition first, and a well-defined EM transition will then allow us to make a 
reliable evaluation of the strong isospin violation mechanism. 

The EM transitions can be studied in the vector meson dominance (VMD) model. Recently, a system- 
atic investigation of the role played by the EM transitions in J/ip(ip') — > VP, where V and P denote light 
nonet vector and pseudoscalar mesons, respectively, was reported in Refs. [H, EH; and the up-to-date 
experimental data provided a good constraint on the VMD model. For <\> — > ojtt , the VMD approach has 
great advantages: on the one hand, the <j> and u> meson masses are very close to the p mass. Hence, the 
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EM form factors can be constrained by the precise data for the p° meson mass and width [2|. On the 
other hand, since other heavier vectors are rather far away from this kinematic region, their contributions 
to the form factor will be limited. The dominant mechanisms can thus be clarified. The availability of 
experimental information for tf> — ► 771"° and pir + 7r + 7r~7r° 14] is also an advantage for quantifying the 
EM contributions. 

The isospin-violating strong decay can be related to the OZI-rule violation at low energies via inter- 
mediate hadronic loops as proposed by Lipkin [l5L [H[ . Microscopic interpretation of such a scenario as 
a mechanism for the OZI-rule violation was investigated by Geiger and Isgur in a quark model [13, 
For instance, an ss pair of 1~ can couple to non-strange nn = (uu + dd)/y/2 via K K, K*K + c.c, etc. 
Suppressions of such an OZI-rule- violating process come from the cancellations between the intermediate 
meson loops and off-shell effects on the intermediate states fl9l I20I] . Qualitatively, at high energies, where 
the mass scale of the intermediate states becomes unimportant, one would expect a "perfect cancellation" 
among all those intermediate states, and it recovers the OZI rule. At low energies, where the mass scale 
of the individual states is dominant, the "perfect cancellation" will break down due to e.g. m u 7^ md 
originated from the chiral symmetry breaking. The OZI-rule violations hence give rise to the recognition 
of isospin symmetry breakings. 

Such a mechanism in (ft — > W7r° decay can be described as follows: In <fi — > unr , the intermediate 
charged and neutral kaon loop transitions are supposed to cancel out if the isospin symmetry is conserved. 
However, due to small mass differences between the u and d quarks, the charged and neutral kaons will 
also have small differences in mass, i.e. m^o —m K ± = 3.972 ± 0.027 MeV [H, and they are coupled to the 
4> meson with slightly different strength. The hadronic loops will then have "imperfect" cancellations and 
lead to measurable isospin violating branching ratios. This drives us to investigate the contributions from 
the intermediate meson exchanges to 4> — > W7r°, which are not only an OZI-rule violating mechanism, but 
also a source of isospin violations. 

A reasonable approach is that at hadronic level, we study the EM and hadronic loop contributions co- 
herently with the aid of the up-to-date experimental data. It will enable us to quantify these two isospin 
violating sources with some obvious advantages: i) At hadronic level, we can extract couplings from 
independent experimental measurements without knowing all the details about the quark distribution 
functions. This technique has been broadly applied to the study of non-perturbative long-range interac- 
tions in the hadronic decays of heavy quarkonia, especially in charmonium decays (2ll. |22|. l23l I24I |25| . ii) 
Adopting the experimental constraints on the meson masses and effective couplings, we also avoid the 
details about how the difference of the u-d quark masses leads to the corrections to the decay constants. 

In the next Section, we first analyze the EM <fr decay in a VMD model and then present our intermediate- 
meson-exchange model with effective Lagrangians. The numerical results for <f> — ► W7r° are given in Section 
III. An extension of this approach to J/ip(tp') — > uitt° is also discussed. A summary is then given in Section 
IV. 



II. THE MODEL 
A. Electromagnetic decay in VMD model 

The Vj* coupling is described by the VMD model [26| , 

C Vl = £ ^V^A" , (1) 

where eMy/fy is a direct photon- vector- meson coupling in Feynman diagram language, and the isospin 
1 and component of the EM field are both included. It should be noted that this form of interaction 
is only an approximation and can have large off-shell effects arising from either off-shell vector meson or 
virtual photon fields. In this approach we consider such effects in the V7P coupling form factor which 
will then be absorbed into the energy-dependent widths of the vector mesons. 
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The typical effective Lagrangian for the VjP coupling is: 

C VlP = 9V \ P f iq2) e^d^d a A^P (2) 

My 

where V v (= p, w, (ft, J/tp, tp' ■ • •) and A 13 are the vector meson and EM field, respectively; My is the 
vector meson mass; e^ap is the anti-symmetric Levi-Civita tensor. The coupling constant gvjp(q 2 ) is 
off-shell and involves a form factor due to the virtuality of the photon. It can be expressed as 

gv^piq 2 ) = g Vl p{^)^{q 2 ) , (3) 

where gy 7 p(0) is the on-shell coupling and can be determined by vector meson radiative decays (T3. 
e.g. uj — > 77T° and (ft — > jir . 

In the VMD model, we can decompose the virtual photon by a sum of vector mesons as shown by 
Fig. [U The amplitude for Process-I (i.e. Fig. Hp)) can be expressed as 

M EM-I _ Mv e gW> a J3 V- v (as 

* ~ L> f v Ml - M v + iM v T v U M, £a ^ p ^ p <P £ 4> < 



v 



where g^v-n is the VVP strong coupling constant, and Ty is the total width of the intermediate vector 
meson. This gives 

gv,p(i 2 ) = gv.p^Ho 2 ) = _ ul\ i Mv T v ' (5) 

which relates the on-shell coupling gv 7 p(0) to an off-shell coupling with form factors. 
Similarly, the transition matrix element for Process-II (Fig. QIII)) can be written as 

M EM-II J MY. e a Prfi u ( ps 

fl ~^fv Ml -M v + iM v r v f w M<j> £ ^P^P* £ * ' W 

where g^y^ is again the strong coupling constant. 

In the VMD framework, it also allows contributions from Process-Ill (Fig. [ljlll)) of which the expres- 
sion is 

M EM-III = \- ^1 M K e e 9ViVi* n a F 0^ F v 

ft M fv 2 Ml - M Vv + iM Vl T Vl Ml - M\ + iM V2 T v , f u U M Vl £ «^ £ ^ £ t> > 

where V\ and V2 are intermediate vector mesons which are different from u> and (ft when they are connected 
to these two states by the virtual photon. However, since we adopt experimental data for (ft — ► p°ir° in 
Process-II to determine the g^po^o coupling, contributions from Process-Ill will have been included in 
Process-II. Nonetheless, we note in advance that exclusive contributions from Process-Ill are negligibly 
small. Therefore, we will only concentrate on the first two processes in this study. 
The following points can be made about (ft — > W7r°: 

i) We argue that the dominant contributions are from p° in this kinematics. Contributions from higher 
states will be relatively suppressed because their masses are larger than the virtuality of the photon. 
Other suppressions from the V^* and VVP couplings are also expected. Basically, those higher vector 
mesons are farther away from the (ft and u> masses than the p°. We thus make an approximation of 
Eqs. and (© by considering only the p meson contributions: 

Mf? = Mf L - 1 +Mf i M ~ n 

= ^e^pZe^e; , (8) 
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where the EM coupling (Jem has a form: 



9 em — -7- 

Jp 



Ml e Ml 

Ml - Ml + l M p T p 9 ^° + /„ Ml - Ml + iM p T p 



(9) 



with r p and r u the total widths of p° and u>, respectively. 

ii) The vector-meson-photon couplings, e/ fv, can be determined by V 



e 



3r, 



-,1/2 



2a e |p e 



where |p e | is the electron three- momentum in the vector meson rest frame, and a e 
fine-structure constant. 



(10) 

1/137 is the 



iii) The coupling, g^ p o^o — 85, can be well determined by either tu — > jtt or uj — > tt e + e~ [2| in the 
same framework. 

iv) For g^pOpa, the KLOE measurement suggests that 4> — > pn — > 7r + ir~ir° has a weight of 0.937 in 
<j) -»■ tt+tt-tt 14]. This gives 



0.937 x Tl 



♦PTT+tt+tt-ttO 12ttM? 



(11) 



with |p| denoting the three- vector momentum of the final state meson in the 0-rest frame. It is reasonable 
to assume g^po^o — g^p+K- — g^p-n+i Thus, the coupling constant can be determined: g^po^o — 0.68. 

On the other hand, the coupling g^pa^o can be extracted in <fi — > jit by assuming that the p° is the 
dominant contribution to the form factor. This leads to 



g<t>p°n a 



\p\ 3 (e/f P ) 2 



1/2 



Ml 



0.68 , 



(12) 



where the p meson width is included. These two results are in excellent agreement with each other 
and highlight the necessity of considering the width effects of the p° pole in the form factor. Also, this 
evidently shows that the p° pole is the dominant contribution in the <f> meson radiative decays, and the 
VMD approach indeed provides a reliable description of the EM transitions in tfi — > lutt . 

In the above treatment all the couplings are determined by experimental data and there is no free 
parameter in the calculation of the EM decay couplings. 



B. Intermediate KK(K*) + c.c. loop 

As discussed in the Introduction that one, in principle, should include all the possible intermediate 
meson exchange loops in the calculation. In reality, the break-down of the local quark-hadron duality 
allows us to pick up the leading contributions as a reasonable approximation (l5l . Il6| . In the <fi meson 
decay, the leading branching ratio is via <p —> KK, which makes the intermediate KK rescattering via 
K* exchange a dominant contribution. Apart from this, tpK*K coupling is sizeable in the SU(3) flavor 
symmetry which also makes the intermediate KK* + c.c. rescattering via kaon and/or K* exchange 
important contributions in 4> —> ljtt . Contributions from higher mass states turn to be suppressed at the 
<f> mass region. We take this as a reasonable approximation in this work, and formulate the contributions 
from i) intermediate KK(K*) loop; ii) intermediate KK*(K) loop; and iii) intermediate KK*(K*) loop. 

The transition amplitude for <f> — > urn via an intermediate meson loop can be expressed as follows: 

M ti = I E ^^hpD ■ (is) 

J (2tt) 4 A-^ aia 2 a 3 

K*pol 
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For KK(K*), the vertex functions are 

igiipi -Ps) ' £0 
^£ a/3 ^pS«£2 (14) 

where 31, 32, and 33 are the coupling constants at the meson interaction vertices (see Fig. [lj. The four 
vectors, p^,, p u , and p^o are the momenta for the initial <fi and final state u> and 7r meson; The four- vector 
momentum, pi, P2, andp3 are for the intermediate mesons, respectively, while a\ = pf — mf, 0,2 — p\ — m\^ 
and 03 = p\ — m\ are the denominators of the propagators of intermediate mesons. 

The form factor T(p 2 ), which takes care of the off-shell effects of the exchanged particles, is usually 
parameterized as 

^ 2 )=(^rzf)"> (15) 

where n — 0, 1, 2 correspond to different treatments of the loop integrals. 

The coupling constants for the charged and neutral meson interactions are denoted by subscription 
"c" and "n" , respectively. In the charged meson exchange loop, coupling g± c can be determined by the 
experimental data for — > K + K~ + c.c, 

2 6lrM 4> 

die = | Pi ^<t>->K+K-+c.c. ) (16) 

where T^ K+K - +C c — (49.2 ± 0.6)% x T to t 2]. For the neutral channel, g\ n is determined by <j> — > 
K°K° + c.c. for which we adopt T 4 ,^ KsKl = (34.0 ± 0.5)% x T tot to derive: 

2 6nM l 

5 i« = Tp — \3 T <f>^ K sK L ■ (17) 

l-r l n | 

The coupling constant g^ c and g^n can be deduced through the decay K* — » Kir. For example, g^ n is 
determined by K*° -> Jf°7r°: 

2 2 67rAf^,o 

53n = 9K*°K°n° = ip|3 r^*°— K°7r • (18) 

It shows that within the precision of the experimental data for K*° — > K°tt° and A"* ± — > if^T, coupling 
9K*°K°Tr° nas the same value as gK* ± K ± ^ ■ The extracted values are listed in Table [IJ 

The relative signs between the couplings are determined by the SU(3) flavor symmetry relations [27| : 

.93c = — <?3ri = 9K*-K-n° = ~9K' + K+-K a = 9K*°K 7r° = -.9x*o/f °7r° ■ (19) 

Note that the above equation is to illustrate the relative signs instead of the values for the coupling 
constants. 

The coupling constant 32 cannot be directly derived from experiment. But it can be related to the 
Lup°ir° coupling via the SU(3) flavor symmetry: 

32c = 92n = 9ujK*-K+ = 9ujK' + K- = 9ljK*"K° = 9ujK*°K° — 3wp 7r°/2, (20) 

where, again, the relative signs between the charged and neutral couplings are determined by Ref. [2?| . 

With the couplings determined as the above, one can see that a relative sign arises between the 
amplitudes for the charged and neutral meson exchange loops. We then distinguish these two amplitudes 
as follows: 

M ti = M% + Mf, , (21) 
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where Mj { and have similar structures except that the couplings and masses involving the interme- 
diate charged and neutral mesons are different due to the isospin symmetry violations. The nonvanishing 
cancellation thus can contribute to the isospin- violating branching ratios. 

To proceed, we treat the loop integral in two different ways. Firstly we apply an on-shell approximation 
(Cutkosky rule) for the intermediate KK, which will reduce the loop integration into an integral over 
the azimuthal angles defined by P3 relative to p^. This approximation picks up the imaginary part 
of the transition amplitude, and with n = 0, 1, 2, we can examine the effects from the form factors. 
Disadvantage of this treatment is that for intermediate mesons of which the mass threshold is above the 
<j> mass, their contributions to the imaginary (absorptive) part vanish though their contributions to the 
real (dispersive) part may be sizeable. Because of this, we also consider the loop integrals including the 
dispersive part in a Feynman integration. To kill the ultraviolet divergences, we include the form factors 
with n = 1 and 2 for a monopolc and dipolc, respectively. Below are the details. 



1. Integrations with on-shell approximation 



By applying the Cutkosky rule to the loop integration, we can reduce the transition amplitude (e.g. 
for the charged meson loop) to be: 



l f i= 32^M-J dn ^i' (22) 
with 



T c ee (T 1 T 2 T 3 ) C = mc9 ^ 93c ^ v eZplM^ ■ P3c . (23) 

The integration is over the azimuthal angles of the momentum p3 C relative to the momentum of the 
final state ir meson. The kinematics are defined as p u = (E u , 0, 0, |P W |), p-n = (E n , 0, 0, — |P W |), and 
Plc = (P3c-Pn) 2 = Ml + m| c -2^3 C + 2|P T ||p 3c |cos0. 

Similarly, we obtain the amplitude for the neutral meson loop: 



M?. = J^L / dfl T f^} , (24) 
fl 32^2 M0 J p\ n - m\ n 

with 



T n = (TiT 2 T 3 )„ = — 4e Q/3 ^£ w p 3 J ,Xp w £ • p 3n . (25) 

Note that the momenta and masses for the intermediate states arc different between the charged and 
neutral cases as denoted by the subscription "c" and "n" , respectively. 

The nonvanishing amplitudes require the vector meson polarizations to be taken as either (e w ,£^) = 
(+, — ) or (—,+). We then obtain 

M fi{ +, ) = -MM-, +) = - 9m9 t^y^ ■ (26) 

where 

Sin2 2 ^?)sin^. (27) 



(i) With no form factor, i.e., T(p\) = 1, the integral becomes: 



2 A 2 -I 1 + A 
^ + -4^ l0g T— 4 



(28) 
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(ii) With a monopolc form factor, i.e., T{p%) = (A 2 — rn^)/(h? — p 2 ,), the integral becomes: 



mi, - A 2 



A 2 -I 



log 



1 + A 



1-B 2 



log 



1 + B 



AB A 2 (A-B) b l-A B 2 (A-B) b 1 - B 
(iii) With a dipole form factor, i.e., ^-(p 2 ,) = [(A 2 — m?,)/(A 2 — P2)) 2 , the integral becomes 



1 = 



(to 2 - A 2 ) 2 
A S B 2 (A-B) 2 



2B(A- B)(B 2 - 1) A 2 -I 1 + A AB 2 - 2B + A l + B 



B 2 (l-B 2 ) 



A 



log 



1-A 



The kinematic functions are defined as 



A. 



M 2 + mf - 2E 1 E l ^ 
B s = M 2 +m 2 -2^K 
A = -2\ Pl \\P u \/A., 
B = -2|pi||P w |/B,. 



ml, 
A 2 , 



B 2 



log 



1 - B 



(29) 



(30) 



(31) 



(32) 



2. Feynman integrations with form, factors 

With the form factors, the ultraviolet divergence in the Feynman integration can be avoided. For the 
charged meson loop as an example, the integral has an expression: 



d 4 p 2c 
(2tt)< 



E 



[igicipic - P3c) ■ £^}[Tr L£ af3^PZ £ uP2c £ 2}[ i g3c(P-K + P3c) ' S 2 



K'pol 

With a monopole form factor, we have 

JVl fi — ^ £ aP l iiyP UJ £ UJ P c f,£ ( l > J ax 

while with a dipole form factor, we have 



(Pic - m \c)(v\c - m 2 3c )(p 2 2c - m 2 c ) 



HpIo) 



A(mi e ,m 3e ,A) 



JVL fi — ~ r £ aPp,yP UJ £ UJ P l p£ ( t > 



y(A 2 - ml c ) 
A(mi c ,TO 3c , A) 

where the function A is defined as 



dx 
^0 



o (Ait) 2 A(mi c ,TO 3c ,m 2c ) 

A(mi c ,TO 3c , A) 



l-x 



dtj 



(An) 2 



log 



A(toi c ,to 3c ,to 2c ) 



(33) 

(34) 

(35) 
(36) 



A(a, 6, c) = Mj(l - x - y) 2 - (M 2 - M 2 - M 2 )(l - x - y)x + M 2 x 2 - (M 2 - a 2 )(l - x - y) 

-(M 2 -b 2 )x + yc 2 . (37) 

Expressions for are essentially the same as with <?i c .2e,3c and m\ c .2c,3c replaced by <7in,2n,3n 
and mi nt 2n,3n, and we do not repeat them here in order to save space. 



C. Intermediate KK*(K) + c.c. loop 
As shown by Fig. [21 the vertex functions for the KK*(K) + c.c. loop are 



= i/i 



l ) 3 t 3 1 



r 3 = i/aCfTT - Pa) • £3 ■ 



(38) 
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where /i,2,3 are the coupling constants and and T(j>\) is the form factor. 

Similar to the previous Section, one finds that a relative sign arises from the charged and neutral meson 



exchange loops, which can be distinguished by Mfi 
charged amplitude with a monopole form factor: 



+ Mf}. Thus, we have the expression for the 



M c _ flchchc a H i> 

Jvl fi — j^j £ aP^Puj e ujP l j > £ 4, 



and with a dipole form factor: 



M c flchchc a fl u „ 

y{A 2 - mfj 



dx 



dx 



dy 



(4tt) ; 



log 



dy 



(4tt) 2 



log 



A(m lc ,m 3c , A) 
A(mi c ,m 3c ,TO 2c ) 

A(mi c ,m 3c , A) 



A(m lc ,m 3c ,m 2c ) 



A(mi c ,m 3c , A) 



(39) 

(40) 
(41) 



In the above two equations the intermediate meson masses 7721,2.3 are from the KK*(K) loops, which are 
different from those in Eqs. (f3~4")> and (|3"5|) . 

In the KK*(K) loop, the coupling constant g$K*K is related to g wp o v o in the SU(3) flavor symmetry: 



flc — fin — g^K' + K 9(f>K*-K+ — 9^K*°K° ~ 9<j>K* a K° ~ 9u)p°-K°/V2 



(42) 



where we neglect the possible differences caused by the isospin violation between the charged and neutral 
channel. The reason is because this loop contributions are negligibly small and such a differences cannot 
produce measurable effects. At the ojKK vertex, the coupling g^Kii can be related to cf>KK by the 
following relation: 



he — g^K+K —9uK-K+ — 9(j>K+K- / V2 , 

hn = 9a)K°K° = ~9uK°K = S^o^o/V^ , 



(43) 



where we assume that the isospin breaking in the luKK couplings is similar to that in the <fiKK ones. 
The absolute values of the coupling constants are listed in Tabic |TJ 



D. Intermediate KK*(K*) + c.c. loop 



We also consider the transition amplitude from the intermediate KK*(K*) + c.c. loop (Fig. [2|, which 
can be expressed the same form as Eq. (IT31 except that the vertex functions change to 



M,:, 

ih-2 



x/3iMvf$ £ $P3 £ 3 > 



Tn = ,ai 1 ,n a F P Pt 1 F v 

J-2 — m2 ta'^'/iVft £ 2 r uj £ uj i 

To = H^-F nan ti „ n a ' ' F^" T)^" F V '' 

J 3 — m3 £ a"l3"n"v"P 2 £ 2 P3 £ 3 



(44) 



where /i,2,3 are the coupling constants and J^ip^) is the form factor. 

Similar to the above Sections, there is a relative sign arise from the charged and neutral meson exchange 
loops, i.e. Mfi = + AF/-, and we only give here the expressions for the charged amplitude with a 
monopole and dipole form factor respectively, 



and 



, , c hxchlchzc a a u „ 

M n = AT m m ea P^P^P4> e 4> 

M < f > m2 C m3 C 



dx 



dy 



l-x-y 2 
dz 



B 



(4tt) 2L Ai 2A? 



2J ' 



M 



fi 



- e a0ni/PZ e uiP4, e d> 



dx 



dy 



l-x-y 



dz 



2 A B 
(4^Af ~ Af J 



(45) 



(46) 
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with 

A = \{2x + ^z~ 1)(M| - Ml - Mlo) + \xMlo + l z M* , 

B = (x + z- l)xz[MlMl a - - Ml - M^f] , 

A 1 = x 2 M 2 + z 2 Ml - xz(M 2 - Ml - Ml,) - z(M% - M\ c ) 

+yM 2 2c - x(M% - Ml) + (l- x - y - z)A 2 . (47) 

In this transition loop the intermediate meson masses mi. 2, 3 correspond to K, K* and k*. Quantities 
^1,2,3 denote the corresponding vertex coupling constants with the relative signs given by: 

hzc = —h>3n = —9k* K*°tv = —9K*°K*°n° — —9k* + K* + tt ( > = —9K*-K—ir° = 9u P °ir°/2 • (48) 

III. NUMERICAL RESULTS 
A. Branching ratios from EM decay transition 

The </> meson EM decay turns to be very sensitive to the p° mass pole and decay width in the VMD 
model. This is because their masses are close to each other. As a test, in the infinitcly-narrow-width 
limit, i.e. T p — — GeV, the branching ratio turns out to be overestimated: BR EM — 1.46 x 10~ 4 , 
which is more than two times of the experimental value. This may not be surprising since one should 
adopt the mass eigenstates in the calculation instead of the isospin eigenstates in degenerate perturbation 
theory. Therefore, we apply the experimental data for the intermediate vector meson masses and widths 
in the calculation. 

With the width of the p meson included, we obtain BR EM = 1.68 x 1CT 5 , with M p = 775.9 MeV and 
T p = 143.9 MeV [H]. With the PDG average, i.e. M p = 775.5 MeV and T p = 149.4 MeV, we have 
BR EAI = 1.67 x 10~ 5 . This explicitly shows an important role played by the p meson. 

We also examine the relative strength between Process-I and II. Their exclusive contributions to the 
branching ratios are BR EM ^ 1 = 1.45 x 10~ 5 and BR EM ~ 11 = 4.56 x 10~ 7 , respectively, which shows 
that Process-I is dominant over II in the (f> decay. 

The above results suggest that the EM transition alone cannot account for the observed branching 
ratio for tfi — * unr°. We hence need to look at the contributions from the intermediate meson exchanges. 

B. Branching ratios from hadronic loop under on-shell approximation 

Under the on-shell approximation only the intermediate KK will contribute since the threshold of any 
other strange meson pairs will be above the <fi mass. 

Without the form factor, the branching ratio from KK(K*) loop is 3.02 x 10~ 6 . This number is 
much smaller than the EM contributions. Apart from the significant cancellations between the charged 
and neutral channel amplitudes, another reason is because of the kinematic suppression on the absorp- 
tive amplitudes, i.e. the intermediate A'A' is close to the <fi mass. Similar phenomena are observed in 
J/tp — > 7/o(1810) — * juxf) at the higher mass tail of the /o(1810) [28|. At least it is reasonable to under- 
stand that contributions from near-threshold intermediate meson rescattering are limited in the on-shell 
approximation. 

In order to investigate the role played by the form factors, we present the calculation results in Fig. [3] 
for three cases: i) The hadronic loop has a dipole form factor (solid curve); ii) The hadronic loop has a 
monopole form factor (dashed curve); and iii) no form factors are included (dot-dashed line). It is easy 
to understand that under the on-shell approximation the calculation without the form factors for the 
hadronic loops will have the largest contributions to the branching ratio. In contrast, the inclusion of a 
monopole form factor suppresses the hadronic loop contributions, and a dipole form factor leads to the 
most suppressions. These three results then converge to the same value when A — > 00 as shown in Fig. [3] 
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The overall results in terms of A including the EM and hadronic loop amplitudes are presented in Fig. 2] 
for two different phases, i.e. on the left panel the EM amplitude is out of phase to the hadronic loop 
(destructive addition), while on the right panel these two amplitudes are in phase (constructive addition). 
On the left panel the horizontal line reflect the largest cancellation between the EM and hadronic loop 
amplitudes with no form factor suppressions. At small A region, the cancellations are small for both 
monopole and dipole calculations since the hadronic loop amplitudes are small in both cases as shown 
by Fig. G3 These three curves smoothly approach the same value at high A where the hadronic loop 
contributions become negligibly small. 

On the right panel the EM amplitude is in phase to the hadronic loop. In the case that no form factor 
introduced in the hadronic loop, the constructive addition of the EM and hadronic loop amplitudes gives 
BR = 2.55 x 10~ 5 . For the monopole and dipole form factor, the constructive effects increase with 
parameter A since the exclusive hadronic loop contributions are small in small A region. It shows by the 
dashed and solid curve that the inclusive branching ratios converge to the dot-dashed curve at large A. 
In this constructive addition, the maximum branching ratio is still smaller than the experimental data, 
which is a sign for the underestimate of the hadronic loop contributions in the on-shell approximation, 
and implies the need for contributions from the dispersive part, i.e. from intermediate mesons above the 
4> mass. 

C. Branching ratios from Feynman integrations 

Note that we are interested in a small effect arising from cancellations between two sizeable amplitudes. 
Since the charged and neutral amplitudes distinguish themselves by the mass differences between the 
charged and neutral particles involved in the loop transition, it makes the behavior of the cancellations 
very sensitive to the choice of the cut-off energies. Again, it is necessary to investigate the A dependence 
of the hadronic loop integrals. We first study the exclusive behaviors of the KK(K*), KK*(K) and 
KK*(K*) loops and then combine them with the EM transitions to study their interferences. 

In Fig. the KK(K*) loop in terms of the cut-off energy A is illustrated. The left panel is for 
a monopole form factor, while the right one is for a dipole type. The dashed and dot-dashed curves 
are contributions from the charged and neutral meson loop, respectively, and the solid curves are their 
differences. In fact, the differences between the dashed and dot-dashed curves are so small that it is 
hard to distinguish them as shown by the figures. Their cancellations leave only a small residue quantity 
accounting for the isospin violation effects. 

The dependence of the details of the cancellations to the cut-off energy turns out to be more dramatic 
with a dipole form factor as shown by the right panel of Fig. [5] Although the integral for both the charged 
and neutral meson loops has a well-defined behavior, details of the cancellations as shown by the solid 
curve has an oscillatory behavior at small A. This is understandable since the difference between the 
charged and neutral meson loop integrals has a complicated dependence on the couplings, and the mass 
differences between the charged and neutral kaon and K* in the propagators. For large A, the integral 
difference smooths out since A becomes the major energy scale. 

In Figj5] there arc dips appearing at small A for both monopole and dipole form factors. This is due 
to the factor A 2 — m 2 Kt in the numerators of the form factors and the largest cancellation between the 
charged and neutral meson loops. 

For the P-wave <fi — > lott decay, the form factor favors a dipole behavior with relatively large A in 
order to account for the off-shell effects. Guided by the solid curve on the right panel of Fig. we 
argue that A ~ 1.5 ~ 2 GeV is appropriate for the hadronic loop contributions. Also, in this region, the 
integral difference has a well-defined smooth behavior. In the case of monopole form factor, to describe 
the experimental data, A must have a relatively smaller value, i.e. < 2 GeV. Otherwise, the branching 
ratio will be overestimated. Due to this ambiguity, we leave the value of A to be determined by the 
experimental data. 

The KK*(K) loop contributions are presented by Fig. [6] for the monopole and dipole form factors. 
Similar to Fig. the intermediate charged and neutral meson loop contributions to the branching ratios 
are compared with each other as denoted by the dashed and dot-dashed curves, while the solid curves are 
given by their amplitude differences. Interestingly, the KK*(K) loop contributions turn out to exhibit 
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a smooth behavior with both monopole and dipolc form factors, and their magnitudes are comparable 
with the KK(K*) loop. Again, the dips are related to the factor A 2 — m 2 K in the numerator of the form 
factors and the largest cancellation between the charged and neutral meson loops. 

In Fig. [3 the A-dependence of the exclusive contributions from the KK*(K*) loop are presented. 
Compared with the other two loops, the exclusive branching ratio decreases in terms of the increasing A. 
As a result, its interferences with other channels around A = 1.5 ~ 2.0 GeV turn to be small. 

Adding the hadronic loops to the EM amplitude coherently, we examine two phases in Fig. [5] in 
terms of the A, i.e. constructive (left panel) and destructive additions (right panel). It shows that with 
A = 1.8 ~ 2.3 GeV, the constructive addition with the dipole form factor for the hadronic loops gives 
the branching ratio in agreement with the experimental data, while with the monopole form factor, A 
requires a range of 1.2 ~ 1.5 GeV. These cut-off energy ranges are consistent with the commonly accepted 
values. For a destructive addition between the EM and hadronic loop amplitudes as shown on the right 
panel, we find that the dipole form factor cannot reproduce the data within A = 1 ~ 2.6 GeV due to the 
significant cancellations between the EM and hadronic loop transitions. In contrast, with a monopole 
form factor for the hadronic loops the destructive addition can still reproduce the data around A = 2.3 
GeV. However, this value of A turns to be out of the commonly accepted range for a monopole cut-off 
energy. In this sense, it shows that the data favor a constructive phase between the EM and hadronic 
loop amplitudes. 

The dipole form factor might be even more preferable. As we have discussed earlier that the P-wave 
decay will generally favor a dipole form factor, we hence argue that the constructive addition between 
the EM and hadronic loop amplitudes with a dipole form factor is a favorable mechanism accounting for 
the experimental observation of BR(<f> — > ujtt ) — (5.2~^\'f) x 10 -5 [2]. In Table HT1 branching ratios of 
the exclusive and coherent (constructively) additions of the EM and hadronic loops with the dipole and 
monopole form factors are listed in comparison with the data. 

In comparison with the results given by the on-shcll approximation, it shows that the dispersive part 
of the loop transitions plays an important role in reproducing the data. 

D. Hadronic loop contributions to the isospin violations in J/ip — > ujtt 

Similar to <fi — > unr°, the decays of J ftp — > ujtt and ip' — > ujtt are also isospin violating processes 
via DOZI transitions. Their branching ratios are measured in experiment, i.e. BR(J/tp — ► ujtt ) — 
(4.5 ± 0.5) x 10~ 4 and BR(ip' -> ujtt ) = (2.1 ± 0.6) x 10" 5 0, which are not significantly suppressed 
compared with J/ip(ip') — > <prj, uirf ', etc. An explanation based on vector meson dominance is provided 
in Refs. [H, [l3[ where the branching ratios are fitted by EM transitions with an appropriate form factor. 
It also shows that Process-I is the dominant contributions to the branching ratio while Process-II is 
negligibly small. In this study, a natural question is about the role played by the hadronic loops and 
their contributions to the branching ratios. 

Interestingly, J/t/j — > K* K is one of the largest decay modes, from which relatively large couplings 
for the J/ipK*K vertex can be derived. However, due to the heavy mass of J/tp, suppressions on 
the loop amplitudes become crucial. With the cancellation between the charged and neutral KK{K*) 
loops, the hadronic loop contributions to the branching ratio turn out to be orders of magnitude smaller 
than the data. In i\j' decay, the cancellation between the charged and neutral KK*(K) loops is not as 
significant as that in J/-0 where the branching ratios, BR(J/4> — > K* + K~ + c.c.) = (5.0 ± 0.4) x 10~ 3 
and BR(J/ip -> K*°K° + c.c.) = (4.2 ± 0.4) x 10" 3 are_ close to each other. In contrast, BR(ip' -> 
K*+K- + c.c.) = (l.Ttg:?) x 10~ 5 and BR{ip' -> K*°K° + c.c.) = (1.09 ± 0.20) x 10" 4 have large 
differences, and have contained significant contributions from the EM transitions (l2.H3|. This favors to 
maximize the isospin violation effects in the hadronic loops. However, due to the suppression from the 
off-shell form factors, the hadronic loop contributions will still be negligibly small compared with the EM 
transitions. 

The numerical calculations show that the branching ratios from the intermediate KK(K*), KK*(K) 
and KK*(K*) loops in J/ip(ip') — > ujtt are orders of magnitude smaller than the data. This result 
suggests that the EM transition is likely the dominant isospin- violating process in the vector charmonium 
decays into light vector and pseudoscalar mesons. Thus, it enhances the argument [H, HH that the long- 
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standing "pn puzzle" in J/ip(tp') — > VP is mainly due to the strong destructive interferences from the 
EM transitions in ip' — ► pir which leads to the abnormally small branching ratio fraction of BR(ip' — ► 
fm)/BR{J/ip -> P tt) ~ 0.2% @. 

IV. SUMMARY 

We investigate the isospin- violating mechanisms in </> — > oj7r° and J/-0 — > W7r° by quantifying the EM 
and strong transitions as different sources of the isospin violations. The EM contribution is constrained 
in the VMD model, and the hadronic loop contributions is studied by relating them to the OZI-rulc- 
violating processes. At hadronic level, the OZI-rulc violations are recognized through the nonvanishing 
cancellations between the charged and neutral intermediate meson exchange loops. In another word, 
the observation of the isospin-violating branching ratios can be viewed as a consequence of coherent 
contributions from the EM transitions and the nonvanishing cancellations among those intermediate 
meson exchanges due to the mass differences between the charged and neutral intermediate mesons and 
different couplings to the initial and final state mesons. 

By extracting the vertex coupling information from independent processes, we can constrain the model 
parameters and make a quantitative assessment of the strong isospin violations via leading KK(K*), 
KK*(K) and A"A"*(A"*) loops. It shows that the dispersive part of the hadronic loop amplitudes have 
important contributions to the isospin violation and they produce crucial interferences with the EM 
transitions though their exclusive contributions are relatively smaller than the EM ones in 4> ~ * W7r ° 
decay. 

We also study the hadronic loop contributions to the isospin violating decay of J/ip(i/j') — > ujtt , and 
find that they are negligibly small. This is consistent with our previous study of the EM transitions in 
J/ip(ip') —> VP, where we argued that the isospin-violating channels, such as ojtt , pr\, pr( and </>7r°, were 
dominated by the EM transitions [H, EH- However, a caution should be given that in J/ip(ip') — > VP 
the s-dependence of the intermediate vector meson widths turns to be a sensitive factor in account of 
contributions from light intermediate vector mesons. A coherent study of e + e~ — > ujtt over a broad 
range of s is thus strongly desired. 
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Coupling constants 


\9</>kk\ 


\9uK*k\ 


\gK*KTr\(\fK*Kn\) 


\ftf>K'K\ 


Charged kaon coupling 


4.49 


4.58 


3.96 


6.48 


Neutral kaon coupling 


4.62 


4.58 


3.96 


6.48 



TABLE I: The absolute values of coupling constants for the vertex interactions. Their relative phases are 
determined by the SU(3) flavor symmetry. 





A (GeV) 


EM transition 


KK(K") 


KK*(K) 


KK*(K*) 


Total 


Exp. 


Dipole 


2.14 


1.66 


0.23 


0.33 


~ 0.0 


5.2 ±0.2 


(5.2+i;i) 


Monopole 


1.38 


1.66 


0.14 


0.56 


~ 0.0 


5.3 ±0.5 


(5.2+^) 



TABLE II: The exclusive and coherent (constructive) contributions of the EM and hadronic loops to the <f> — > loty 
branching ratios with a dipole and monopole form factor. The experimental data is the world average given by 
PDG2006 Q. The branching ratios in columns 3-8 have a unit of 10~ 5 . The errors estimated in column 7 arc 
due to the precisions taken for the exclusive branching ratios. 
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(I) (II) 
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(III) 



FIG. 1: Schematic diagrams for the EM transitions in <j> — » oj7r°. 




(a) (b) (c) (d) 

FIG. 2: Schematic picture for the decay of <^> — > tJ7r via KK{K), KK*(K) and KK*(K*) intermediate meson 
loops. 
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FIG. 3: The A-dependence of the KK(K*) loop contributions in the on-shell approximation. The dot-dashed, 
dashed and solid curve denote different considerations for the form factors, i.e. no form factor, monopole and 
dipole, respectively. 




A (GeV) A (GeV) 



FIG. 4: The A-dependence of the sum of the EM and KK{K*) loop amplitudes in the on-shell approximation. 
The left panel indicates results for a destructive addition and the right panel for a constructive addition. The 
solid, dashed and dot-dashed curves denote different considerations for the form factors, i.e. dipole, monopole 
and no form factor, respectively. 
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FIG. 5: The A-dependence of the KK(K*) loop contributions to the branching ratio in the Feynman integration. 
The left panel indicates results with a monopole form factor, and the right one with a dipole form factor. The 
dashed and dot-dashed curves are contributions from only charged and neutral meson loop, respectively, while 
the solid curves are the results after cancellations between the charged and neutral amplitudes. We note that the 
dashed and dot-dashed curves are close to each other and difficult to distinguish them by sight. 




FIG. 6: The A-dependence of the KK*(K) loop contributions to the branching ratio in the Feynman integration. 
The notations are similar to Fig. [5] Again, we note that the dashed and dot-dashed curves are difficult to 
distinguish by sight. 
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A (GeV) A (QeV) 



FIG. 7: The A-dependence of the KK* (K*) loop contributions to the branching ratio in the Feynman integration. 
The notations are similar to Fig. [5] Again, we note that the dashed and dot-dashed curves are difficult to 
distinguish by sight. 




FIG. 8: The A-dependence of the constructive (left panel) and destructive additions (right panel) between the EM 
and hadronic loops. The dashed curves denote the results for adopting a monopole form factor for the hadronic 
loops, while the solid curves for adopting a dipole form factor. 



